Universality of Regge and vibrational trajectories in a semiclassical model 
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The orbital and radial excitations of light-light mesons are studied in the framework of the domi- 
nantly orbital state description. The equation of motion is characterized by a relativistic kinematics 
supplemented by the usual funnel potential with a mixed scalar and vector confinement. The in- 
fluence of finite quark masses and potential parameters on Regge and vibrational trajectories is 
discussed. The case of heavy-light mesons is also presented. 

o 

O . I. INTRODUCTION 

(N 

^ • One of the most striking feature of the light-light meson spectra is certainly the existence of the so-called "Regge 

Q trajectories" , that is to say the linear dependence of the square mass as a function of the total spin quantum number. 
One can see in Fig. IIVI the remarkable linear trajectories of the p, K* and (j) families. It is well known that this 
behavior can be obtained in potential models by using a confinement potential well adjusted for the kinematics. 
In non-relativistic descriptions of meson spectra, a r^/"^ confining potential must be used 0, while in relativistic 
kinematics the linear Regge trajectories are obtained with a linear confinement 0] and results also naturally from the 
relativistic fiux tube model of mesons Q . It is worth noting that this mass behavior is also expected for heavy-light 
mesons. 

An interesting theoretical approach of light-light and heavy-light mesons is given by the dominantly orbital state 
(DOS) description |5j| for which the leading Regge trajectory is a "classical" result while radially excited states can 
P\J ' be treated semiclassically. In a recent paper 0, this DOS model has been applied to study light-light and heavy-light 
T-H meson spectra with a confinement potential being a mixture of scalar and vector components. In the case of one or 
\l ' two quarks with vanishing masses, it was shown that linear Regge (orbital) and "vibrational" (radial) trajectories are 
, obtained for an arbitrary scalar-vector mixture but that the ratio of radial to orbital energies is strongly dependent 
' on the mixture. 

Oh! In this paper, we extend the result of Ref. 6] for light-light and heavy-light mesons by considering the presence of 
^ ^ quarks with finite masses and the introduction of a Coulomb-like interaction and constant potential in addition to 
I the confinement. In the limit of small masses and small strength of the short range potential, or in the limit of large 
1 angular momentum, all calculations can be worked out analytically. One could expect that the small quark masses 
■ and Coulomb-like potential do not need to be taken into account to study high orbitally excited state. But we show 
here that such effects are interesting to point out. The square mass formula for light-light mesons composed of two 
identical quarks as a function of quantum numbers and potential parameters is established and discussed in Sec. ITU 
J-j ' A comparison of our theoretical predictions with data is performed in Sec. lIIII Some concluding remarks are given in 
^- • Sec.HVI 



II. THE MODEL 



A. basic ideas 

The starting point of our model is essentially the same as the one presented by Olsson Consequently, we just 
remind the basic ideas and skip a lot of technical details which can be found in that reference. In all our formulae, 
we use the natural units h = c — 1. 

We consider a system composed of two spinless particles with the same mass m, which are allowed to rotate 
around their center of mass fixed at the origin, and to vibrate in the radial dimension. The Lagrange coordinates 
are the interdistance r and the rotation angle ip, angle between an arbitrary fixed direction and the line joining the 
two particles. The free relativistic Lagrangian is very easy to write and depends on the Lagrange coordinates and 
their time derivatives. In order to handle the dynamics between the particles, we introduce phenomenologically an 
interaction potential which depends on the radial coordinate r only. In fact, there exist two Lorentz structures for 
this potential: a scalar term S{r) which simply modifies the mass (to — )■ to + S{r)/2), and a vector term V{r) which 
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is simply subtracted from the Lagrangian. Consequently the starting Lagrangian is given by 

L{r,(p,r,(p) 



i2m + Sir)) Wl 



1 



1 



Vir). 



(1) 



This formula can be applied to particles with spin in a very crude manner by considering spin dependent potential 
terms. One can see that L is independent of ip. This means that the angular momentum J = ^ is a good quantum 
number. Taking into account this condition, the Hamiltonian of the system can be written as 



H{r,pr, J) = d4p^ + ^ + (2m + S{r)f + V{r). 



(2) 



The idea of the model is to make a classical approximation by considering uniquely the classical circular orbits (lowest 
energy states with given J), defined by r=rQ, and thus f = and Pr—0. Let us denote 



Mo(J) = i?(r = ro,p,. = 0, J). 



(3) 



By using the condition ^{r ~ ro) ~ and performing a number of analytical transformations (see Ref. Q), one finds 
the fundamental equations (with the conventions 5*0 = S'(ro), and so on) 



Mo (J) = 



ro 



,2 , , (2m + gp) 



(2m + 5o) + S'^ro] 



(4) 



and 



( K,' + \ Vi' + 4 ^^"^ + ^°^ [(2m + So) + S'.ro] 



{2m + So)S'o 



, 1/2 



(5) 



Equation ^ (which is in general a transcendental equation) allows to calculate ro{J). Once this value is obtained, 
it is put into Eq. to get Mq{J). Up to now, we are in position to obtain only the ground state energy for a given 
J. In order to get the radial excitations, it is useful to make a harmonic approximation around the previous classical 
orbits. Here again, the details can be found in Ref. The harmonic quantum energy is given by 



0(J) = 2- 



12 P 

(Mo( J) - Vo) Vo" + + (2m + So)S'^ + 



^0 



1/2 



Mo (J) - Vo 

Then the mass of the system with an orbital excitation J and a radial excitation n (0, 1, 
approximation, by 

M( J, n) = Mo (J) + n{J) (n + i 



(6) 



) is given, within this 



(7) 



In this general formulation, one sees that there exists, in principle, a coupling between the orbital and the radial 
excitations. In relativistic descriptions of two particle systems (for example in the meson sector), it is generally the 
square of the mass which is considered as the function of the orbital and radial quantum numbers. If one assumes 
that ri(J) ^ Mo( J) (see Ref. (S]), then one can write 



M2( J, n) = A/q ( J) + MoiJ)niJ){2n + 1). 



(8) 



B. The meson case 



1. The potential 



In the previous subsection, we gave a formalism applicable to any type of potential. Such a study could be 
undertaken from the numerical point of view, by solving the equations giving tq and reporting the value in the 
expression of Mq and Q. However there are some drawbacks to proceed that way. First, because of non linearity. 
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there can exist several solutions for Eq. lO and it may be difficult to discriminate between them. Second, the physical 
analysis is less obvious when the results originate from a numerical treatment. So, it is highly desirable to study a 
model for which an analytical solution is available. 

In Ref. the author studied the spectra of the mesons or, more precisely, the Regge trajectories for a system 
composed of a quark and an antiquark with to = interacting via a linear confinement potential that is partly scalar 
and partly vector. He showed that a linear behavior Mq (J) cx J automatically emerges. 

In this paper, we revisit this study, but in a more general framework. First, we relax the constraint of zero mass and 
consider the mass of the quarks m as a free parameter. Second, we admit a long range linear confinement supplemented 
by a short range Coulomb part and constant potential. Indeed quantum chromodynamics tells us that it should be 
so. The one-gluon exchange process gives rise to a Coulomb term of vector type, while multigluon exchanges are 
responsible for the linear confinement. Since we ignore the Lorentz structure of the confinement, we suppose, as in 
Ref. that it is partly scalar and partly vector. The importance of each one is reflected through a parameter / 
whose value is for a pure vector, and 1 for a pure scalar. 

Thus, in this paper, we apply the previous formalism, with potentials given by 

S{r)^far (9) 
in which a is the usual string tension, whose value should be around 0.2 GeV^, and 

V{r) = {1- f)ar - - (10) 
r 

in which k is proportional to the strong coupling constant as- A reasonable value of k should be in the range 0.1 to 
0.6. We will show below that we are able to obtain an analytical solution for this rather realistic situation. 



2. Scaling properties 

In our case, we are faced with a number of parameters, whose origin is very different 

• the mass m (in GeV), the string constant a (in GeV^) and the Coulomb constant n (dimensionless) are dynamical 
parameters; 

• the / parameter (dimensionless) reflects our ignorance about the Lorentz structure of the confinement; 

• the quantum numbers J and n of the state. 

In what follows, we will be interested by the variation of M'^{m,a, f, J,n) as a function of quantum numbers J 
and n, and for a given set of physical parameters m, a, k, /. Mq does not depend on n, as well as $7, so that the n 
dependence is rather obvious (see Eq. ||SJ)). Let us focus on the dependence on J and on the physical parameters. 
Actually, it is possible to exploit a scaling law in our equation. One can remark that ^/a is an energy unit, and 
a length unit, so that it is advantageous to work with these natural units. Let us define reduced dimensionless 
quantities, labeled with an overhead bar, relative to these units. Explicitly, we have 

_ m — Mq — M - n _ I- 

m = — ; Mq = —j= ; M ^ — ; II ^ — ; ro = roVa. (11) 
y/a y/a y/a yja 

If we put these values into our fundamental Eq. jSJ, with application to our particular potentials (|9I10|I . we find, after 
some manipulations, the equation for tq 

/2(2/ - 1) + 4m/(2/ - 1) + 2 {2ln\2! - 1) - Kf{\ - /)] rl - 8Krnf{l - /) 
- 4(1 - /)((1 - /)./2 + 2Km2) + f(^K^ + 8J2) - Amf {k^ + 4J^) rl 

-4k [kto^ + 2J2(1 - /)] rl + 4J2 (4j2 - k^) = 0. (12) 

Equation I^J giving the value of Mq leads in this case to 

— (1 - 2/) - 2mfrf, + 4 - 
ro {{^ - f) rl + k) 
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Similarly, Eq. © giving the value of leads to 

2 { [(1 - /)(2/ + MV - 1) ^ + 4Atm/ rg + 3(1 - f){U^ - k") + (4j2 - ^)k] 

[(1 - f)rl + n]}'^' I {r, {-fr% 2lnfrl + U'] } . (14) 
Some comments are in order. 

• We see that the use of this scaling law allows to reduce the number of parameters to four (m, k, f, J) in the 
dimensionless quantities. 

• If both m=0 and k—0, it remains only the powers (0,4,8) in Eq. H12|l so that the equation is always soluble by 
using Tq as a variable; one can see that only one relevant physical solution exists. It should be stressed that this 
case corresponds to the study of (Ref. [^); we checked that our general formalism applied to that peculiar case 
gives the solution presented in that reference. 

• The complicated equation (jsj simplifies in our particular case to a polynomial of eighth order in tq (Eq. idJ). 
In principle, there could exist 8 solutions, but we must reject non physical complex and real negative values, 
and the real positive values which lead to negative values for Mq and Q. Starting from the unique physical 
solution ro(rn = 0, k = 0), we can expect, by continuity, that there still remains only one physical solution for 
small values of m and k. Although we are not able to prove that the physical solution is always unique, it seems 
that this is a general property; we remarked that this was always true for a wide range of parameters and for 
the cases for which an analytical solution can be found (see below). 

• The linear term in f o is always missing in all expressions. 

• The special values f — (pure vector confinement) and / = 1/2 (half scalar-half vector) lead for Eq. (|12|) to a 
polynomial of 6th order only. Moreover, in the case f ~ 0, the polynomial is of 3th order in Tq. 

• Odd powers of Tq disappear in Eq. (|12(l when m — 0. By choosing rg as a variable, we are left with a 4th order 
polynomial. 

Although the things look quite nice now, it is possible to simplify further, using another scaling law. As one will see, 
this new scaling law is very powerful and allows to derive universal conclusions concerning Regge trajectories. Indeed, 
it is possible to eliminate the J dependence with help of the following changes of variables 

X :^ro J-^/^ ; u = mJ-^^'^ : V = kJ-\ (15) 

Using these variables, our polynomial (|12|l is rewritten 

f{2f - 1) x« + iuf{2f -l)x' + 2 [2u\2f - 1) - v,f{l - /)] a;^ - 8vuf{l - /) x'' 
- [4(3/2 - 2/ + 1) + 8(1 - f)u\^ + fv^] x'^ - 4uf{v^ + 4) x^ 

-Av [vu^ + 2(1 - /)] x^ + 4(4 - v^) = 0. (16) 

The interest of this equation is that the variable x now depends only on 3 variables x = x[f, u, v). Alternatively, one 
can put Eq. 113|l under the form 



Mo = J^/^gif,u,v), 

x[{l - f) x'^+v] 



(l-2f)x^-2ufx^ + A-v^ 
g{f,u,v) = „ , , (17) 



and equation (|14l) under the form 

n^2j-^/^h(f,u,v), 

{ [(1 - /)(2/ ~l)x^ + 3v{2f -l)x^ + Avuf x^ + 3(1 - /)(4 -- v^) + (4 - v^)v] 

[il-f)a:' + v]y/'/{x [-fx*-2ufx'+4]}. (18) 

Under these forms, the formalism is the most general and the simplest possible. Except very specific cases that we 
have studied, it is not possible to get an analytical solution for this general formulation. To proceed further, we need 
to make some kind of approximations. 
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e and vibrational trajectories 

From now on, we are interested only in Regge trajectories which give the behavior of in term of J for large value 
of J. It is thus natural to consider the values of u and v defined by Eq. H15|l as small quantities and to develop g and 
h functions, given by Eqs. ((T7|l and ((TH|) . up to order that is to say keeping constant terms, terms with u, terms 
with and terms with v. Such a calculation is very painful and we have been helped a lot by the MATHEMATICA 
software. We skip below all the cumbersome details to retain only the important things. Basically we expand x as 

x^a{f)+mu + l{f)u^ + S{f)v. (19) 

We put this expression into Eq. up to the given order and identify to zero to get the expressions of /?(/), 
7(/) and S{f). Then we report the resulting value of x in g(f,u,v) and h(f,u,v), limited again to the desired order 
to get Mq and il. Lastly we used these quantities in Eq. to obtain M. It is important to point out that the 
expression (|SJ) is valid up to our order (the term in Sl^ is of higher order). At the end of this study, the behavior for 
large J can be traced back, after reintroducing dimensioned quantities, under the exact form 



M 2 = a A{f) J + B{J) m Vol + C{f) + a D{f)K + a E{f){2n + 1) + 0( J-^/^). (20) 

2g(/) 



The coefhcients A, B, C, D and E can be obtained analytically. Another coefficient R{f) — ^^rfr is specially 



important and will be discussed below. Let us define the auxiliary hmctions s{f),t{f)^y{f) by 

s{f) - (1^/^^/2 = 1- 2/ + 3/2, 

Kf) - V<f) + QP, 

vUf = , (21) 



Our coefficients are given by 



s(/) + (l-/)i(/)- 



A{f) - ^[t + 3(l-/)f , 

B{f) - |[(l + /)(3/-l)+t(l-/)], 

C{f) - -^[t(s + /2) + (l-/)(2/-l)]. 



D{f) = -[< + 3(l-/)] 
E{f) - A{f). 



We plot the various coefficients in Fig. IIVI All these quantities are monotonic functions of / and their ranges from 
/ = to / = 1 are given below 

• 8 > A{f) > 4; 

• < B{f) < AV2; 

• 8 > C{f) > 3; 

• -4 < D{f) < -2^2; 

• 4^/2 > E{f) > 4; 

• V2<R{f ) < 2. 

Equation (|20|l is the exact expression valid for large J; we succeeded in getting an analytical expression in our model. 
It is very powerful since it is universal and the various dependencies are explicitly seen. One can make the following 
comments. 
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• The dominant term for large J is linear, so that we recover pure Regge trajectories. As far as the string tension 
is consider as a constant, the slope of the Regge trajectories aA{f) is universal in the sense that it is independent 
of the system (it is also independent on the strong coupling constant). Only the confinement drives the behavior; 
moreover the slope depends on /, the percentage of scalar confinement. 

• The slope of the vibrational trajectories 2aE{f) presents the same characteristics. Even more important, the 
ratio of the slopes for radial and orbital trajectories R{f) depends only on /. This ratio is greater than 1, as 
indicated by experimental data. 

• The next to dominant term behaves like \/j and deforms the straight trajectories for low values of J. It must 
be emphasized that this term is due only to a non zero mass and is independent of the strong coupling constant. 
Moreover, this term is absent whatever the system if the confinement is purely of vector type, since -B(O) = 0. 

• The displacement of the trajectories from the origin has three contributions; one C(f)m? is due to the mass, 
another aD{f)K is due to the strong coupling constant and the last one aE{f) reflects the zero point motion of 
the harmonic vibration. The position above or below the origin depends upon the relative importance of each 
contribution (do not forget that C and E are always > but that D is always <0). The zero point energy of 
the orbital motion cannot obviously be calculated in our model. 

• There is no coupling between orbital and radial motion for large J values (absence of terms nJ). This is only 
a consequence of the Coulomb+linear nature of the quark-antiquark potential. This may not be true for other 
types of potentials. 



It is well known that, in traditional spectroscopy relying upon Schrodinger or spinless-Salpeter equation, it is 
necessary to add a constant term to the potential in order to get the absolute values of the spectra. In those models, 
the effect of this constant is just to shift the absolute spectrum keeping the same relative spectrum. One can raise 
the question of adding a constant potential in our model. We will see that it is not so difficult to answer this question 
in our framework. 

In principle it is possible to add a constant Cy to the vector potential and a constant Cs to the scalar potential. 
Let us discuss separately these two cases. 

If we add a constant Cy to V{r), this does not affect the functions V'{r)^ S{r), and S'{r). From Eq. (|SJ), it is clear 
that the value vq remains unchanged, and thus all the quantities depending on tq remain unchanged except Vq which 
is change to Vo + Cy. Consequently, from Eq. Q), Mq is changed to Afo + Cy and fl remains unchanged since it 
depends only on Mq — Vq (see Eq. ©). Then M is changed to M + Cy. Squaring this quantity and keeping only the 
terms to good order, the net effect of Cy is to add supplementary terms to the B and C coefficients. 

Let us now add a constant term Cs to S{r); V{r) and all the derivatives are unchanged. Since, in Eq. Q, J 
depends on S{r), ro does depend on Cs. But we see that it is always the value (2m + S{r)) which appears. If we 
introduce a new mass m' = m + Cs/2, we see that the new value of rg is simply rolm' , J) instead of ro(m, J). Making 
this replacement in Mq and O, we see that we have now Moirn' , J) and Q{m' , J) instead of the non-primed values. 
So, the only modification in our formalism is just a change from m to m' . But in doing so, we must consider the new 
value of u defined with the modified mass m' for applying a limited expansion in term of J (m! j ^fa can be small with 
respect to 1, while m can be large with respect to y^). Assuming that such an expansion is justified, the only effect 
of Cs is to modify the terms in the ^fj and the constant contributions. 

One can now gather both effects into a single formula. Let us define the dimensionless quantities £ = 1 + |^ and 
77 = The only effect of adding the constants Cy and Cs is to change, in the universal equation (j^OJ), the -B(/) 
and C(/) coefficients by new ones -B'(/, e, 77) and C'(/, e, 77) defined by 



C. 



Addition of a constant term 



i?'(/,e,77) = ei?(/) + 277^1(7), 
C'(/,e,77) = C(/)+ + 26 77- 



^-(1-/) 
2/ 



(23) 



In particular, the A, D and E coefficients are unchanged, so that the slopes of the orbital and radial trajectories 
are still unchanged in this more elaborated potential, and appears really as universal quantities depending only on 
scalar- vector mixture in the confinement. 
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D. Heavy-light mesons 

The same formalism can be applied to the case of heavy-light mesons in the limit of infinite heavy quark mass, as 
mentioned in Ref. 0. In light-light mesons considered here, two identical quarks orbit at the same distance r/2 from 
the center of mass. For a heavy-light meson, the light quark orbits around the fixed infinitely massive quark at the 
distance r. To obtain the heavy-light Hamiltonian, one must perform the following replacements in the light-light 
Hamiltonian (j^J 

r 

• 2m + S{r) -^171 + S{r); 

• V{r) -> V{r). 

After calculation, one obtains 

= ^ A(f) J + ^ B{f) mV^+ C{f) ^ + aD{f)K + a E{f){2n + 1) + 0( J-^/^). (24) 

III. DISCUSSION OF THE MODEL 

Formula H2U|I gives the meson square mass dependence as a function of parameters. One can ask whether it is 
possible to extract from this relation values for the parameters by comparing with available data. The value generally 
considered for the string tension a is around 0.2 ± 0.03 GeV^. Taking to„ « rud « 0.200 GeV and TOs « 0.450 GeV, 
we can see that the parameter m runs from 0.5 to 1. The Coulomb strength k is not exactly known, but in potential 
models a range from 0.1 to 0.6 is currently accepted. So we have to check that formula 120|l can be applied for this 

range of parameters. In Fig. IIVI we compare the reduced square mass calculated exactly by formulae H8I12I14() 
with the one given by Eq. (|20|l . as a function of the total angular momentum J for some values of /, m, k and with 
Cs = C'v = 0. One can expect that the quality of the approximation is better for small rn and k parameters. This is 
precisely what we get. We present the figure for m = 1, a value already large and, nevertheless, the approximation 
still works reasonably well. For / = 0, the approximate value is a pure straight line (-B(O) = 0). The deviation is 
maximum for large n values, as expected. For / = 1, there is a small deformation both for the exact solution and 
the approximate one. Moreover, the curvature is the same in both cases. Note that the concavity of the curve may 
be inverted by introducing non- vanishing values for Cs and/or Cy. For large values of J, it is apparent that our 
approximate expression is always good. For small values of J, the deviation is maximum for k = 0, but we are in 
a region (m large and J small) for which we expect some error. By continuity, other values of / give intermediate 
situations. 

We see that, even in the rather unfavorable case, the approximate expression does not differ so much from the exact 
results. Consequently, we will base our further discussion on Eq. (|20|l . This last equation can rewritten in the form 

M'^{J,n)=AJ + B^ + C + Vn (25) 

with A, B, C and V depending on parameters /, a, to, k, Cs and Cy. Knowing three orbital excitations, for instance 
M'^{J, n), M'^{J + 1, n) and M'^{J + 2,n), it is possible to calculate A and B quantities. Adding a radial excitation, 
for instance M'^{J,n + 1), the quantities C and T> can also be calculated. The formulae are given by 

A = Ki,o{J)M^{J + 2,n)~K2.o{J)M^iJ + '^,n) + K2.i{J)M\j,n), 

B = A{J)-^ [2M^{J + l,n) - M^{J + 2,n) - M^{J,n)] , 

V = M^{J,n + l) - M^{J,n), 

C = KifliJ) y/j(J + 1) M^{J + 2, n) + K2,i{J) ViJ + '^)iJ + 2) M\J, n) 

- K2,o{J) VJ{J + 2) M^{J + 1, n) - n [m\J', n' + 1) - M^{J', n')] , (26) 

where 

A(J) = 2V7TT - V7T2 - VJ and K^j{J) = ^'^ ITt^'^ ~- (^7) 
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The interest of such formulae is that they are universal in the sense that the quantities A, B, C and "D are independent 
on J and n. Given a flavor sector, they could be, in principle, checked for various orbital and radial multiplets. 
Moreover, the quantities A and V are even independent on m, so that they can be checked independently in several 
flavor sectors. Finally, the ratio V/A depends on / only and thus can provide a strong test. 

In principle, one could obtain physical quantities from the experimental masses. Starting from V/A — R{f), one 
can determine the value of / (see Fig. lIV(l . The quantities A{f), B{f), C{f), D{f) and E{f) can then be calculated. 
The A and V quantities can provide us with cross-checked values of the string tension. If Cs = Cy = 0, the B term 
gives the quark mass m; if Cs 7^ and Cy 7^ 0, only the quark mass difference in two flavor sectors can be obtained. 
Calculating parameter k is more involved. 

Unfortunately, our model is only valid for high values of J, and, for such quantum numbers, a few ground states 
arc known and none radial excitation 0. Moreover, uncertainties exist about meson masses. Nevertheless, we can 
try to obtain an estimation of parameter / by considering the mesons of the p and </> families (see Fig. IIV|) . 

The meson masses used to perform calculations are center of gravity of meson multiplets whose members are 
characterized by an internal total spin quantum number equal to one. They are calculated, as well as the corresponding 
uncertainties, with the procedure given in Ref.fs'l. Using the only three available masses J = 1, 2 and 3 for the p family 
and taking into account the error on these data, we can determine the range possible for the function A — aA{f). 
Assuming usual values for the string tension a, we can deduce the possible range for /. Unfortunately, we obtain no 
constraint on this parameter since values from to 1 are compatible with data. The same calculation done for the 
family leads to the conclusion that / must be very close to 0, that is to say that the conflnement is only of vector-type. 
We have remarked that the results are very sensitive to the value of the masses. If other states with higher values of 
J were known, our conclusion could be changed. 

It is also possible to calculate / by considering a meson and its radial excitations {V — 2aE{f)). This procedure 
is more questionable since only the first radial excitation of the J — I meson is known in the p and families. 
Nevertheless, if the calculation is performed for the two families, we find that all values for / from to 1 are 
compatible with the data. 



IV. CONCLUDING REMARKS 



We have shown that light-light and heavy-light mesons exhibit linear orbital and radial trajectories in the dominantly 
orbital states (DOS) model. Slopes of both type of trajectories depend only on the string tension and on the vector- 
scalar mixture in the confinement potential. From our work, it turns out that small finite quark masses do not 
alter significantly the linearity, specially in the case of dominant vector-type confining potential. As expected, the 
Coulomb-like potential has no effect on trajectories but its influence on the meson masses has been determined in the 
approximation of the DOS description. Lastly, we have shown that the ratio of radial to orbital energies depends on 
the scalar- vector confinement mixture only. 

Formula (|20|l gives the parameter dependence of square mass in the limit of small masses and small strength of the 
short range potential, or in the limit of large angular momentum. Outside these limits, we have remarked that this 
formula only differ slightly from exact results in a large range of parameters. This allows to apply our approximation 
for physical situations and to compare our calculations with experiment. 

From data available, it is difficult to determine the value of the parameter /, that is to say the scalar-vector mixture 
in the confinement. Strictly speaking, our results are only semiclassical ones, and experimental radially excited meson 
masses are only known for small value of total spin. Nevertheless, our calculations favor a dominant vector-type 
confining interaction. 

In our study, we completely neglect the quark spin. Actually, this is not an important drawback. Regge and 
vibrational trajectories concern only mesons with an assumed internal total spin quantum number equal to one 
In usual models the spin-dependent part of the potential is the hyperfine interaction stemming from the one-gluon 
exchange interaction (and may be from the vector part of the confinement) In other models, this spin-dependent 
part stems from an instanton induced interactions 0, |^. In both cases, the contribution of the spin to the meson 
masses is small or vanishing for mesons with total spin equal to one. 

In Sec. Ill Dl we mentioned a method to study the heavy-light mesons in the limit of an infinitely massive heavy 
quark. Such a method is used in Ref. 0. It could also be interesting to develop the calculations in the more general 
case of two finite different masses. 
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FIG. 1: Square mass of some light-light mesons as a function of their assumed total orbital angular momentum The 
isovector nn mesons (p family) are indicated by square boxes, strange mesons {K* family) by triangles and ss mesons {(j) family) 
by circles. The mesons represented are characterized by an internal total spin quantum number equal to one. Excepted for the 
I — states, each mass is given by the center of gravity of mesons with the same I but with total spin different (see Sec. lilH . 
The remarkable linear dependence is pointed out by straight dashed lines. 
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FIG. 2: Coefficients giving the square mass I2()|l as a function of parameter /. The ratio R{f) = 2E{f)/A{f) is also given. 
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FIG. 3: Comparison of the reduced square mass M given by formulae 18112114(1 (exact), and the approximate one given by 
Eq. 12UII (app.), as a function of the total angular momentum J for values of /, m and k equal to and 1, and Cs = Cv = 0. 



